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Abstract: This work is the result of the ideas developed by Ken Yoshida about the 
possibility of extending the range of applications of the matrix model approach to the 
computation of the holomorphic superpotential of the /3-deformed = 4 super Yang-Mills 
theory both in the presence of a mass term and in the massless limit. Our formulae, while 
agreeing with all the existing results we can compare with, are valid also in the case of 
spontaneously broken gauge symmetry. 
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1. Introduction 

Understanding the precise infrared behaviour of asymptotically free gauge theories is a 
daunting task. For instance, determination of the confining properties of the theory has 
proven to be impossibly difficult. The reason for the difficulty in analysing the low-energy 
dynamics of the theory is related to the fact that perturbative methods have a limited range 
of applicability as the coupling is growing larger in the infrared. On the other hand, we lack 
systematic non-perturbative methods. Up to date, the most powerful tool at our disposal 
is undoubtedly represented by lattice simulations. Still this is a numerical approach and 
any possible piece of analytic information would be greatly valued. 

Interestingly, however, there exist relevant exceptions to the previous picture repre- 
sented by super symmetric theories. The constraints on the structure of the theory imposed 
by supersymmetry are so strong that it is possible to obtain detailed low energy information 
about holomorphic quantities, such as the holomorphic superpotential. Holomorphicity, in 
fact, protects physical quantities from quantum corrections and/or allows to safely control 
them providing numerous and very useful non-renormalization theorems |||, 

In recent years, there has been a renewed interest in the field after Dijkgraaf and 
Vafa conjectured that holomorphic quantities in the low-energy regime are related to the 
free-energy of an hermitian Matrix Model (MM). The correspondence has passed various 
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non-trivial tests |^ and it has become a powerful tool for studying the low-energy limit 
of A/" = 1 gauge theories with massive chiral superfields. 

At the same time, a particular M = 1 model, known as the /3-deformation of the 
A/" = 4 Super Yang- Mills (/3d SYM), attracted the attention of the community. Leigh and 
Strassler discovered this model in the nineties Q and various aspects of the model have 
been studied ^, but only after its gravity dual was found [l^ the /3d SYM was 
extensively investigated in numerous papers [0, 12, |13|, [l^, |15|, p^ ]. 

In this paper we apply the Dijkgraaf-Vafa correspondence to the study of the holo- 
morphic low-energy superpotential of the /3d SYM model. In particular, we are able to 
compute the low-energy coupling constants controlling the surviving massless degrees of 
freedom, after the \]{Nc) i— s- JliLi U(-/Vi) spontaneous symmetry-breaking has occurred. 
Along the way, as a side result, we also compute the low-energy effective superpotential of 
the massive version of the (3d SYM model as an inverse mass power expansion. 

The paper is organized as follows. In section |2| we review the standard lore about 
the MM approach and the general properties of various deformations of the AA = 4 SYM 
model. In section ^, relying on MM techniques, we compute the holomorphic low energy 
superpotential in various situations. In section ^ we investigate the phenomenon of spon- 
taneous breaking of the gauge symmetry. We end in section ^ with some conclusion and a 
few considerations about possible lines of future developments. A few more technical issues 
are discussed in appendices. 



2. Preliminary material 
2.1 The Matrix Model 



It has been suggested [17| that for a wide class of supersymmetric gauge field theories there 
exists a deep connection between the "free-energy" of certain (zero-dimensional) MM's and 
the holomorphic superpotential as function of the glueball superfield degrees of freedom. 

In its original formulation it was believed that such a correspondence was limited to 
the perturbative corrections to the low-energy superpotential, thus explicitly excluding the 
celebrated Af = 1 Veneziano-Yankielowicz (VY) superpotential ||l^, ^ |^], Wyy- This 
deficiency was attributed to the fact that the overall constant of the matrix model integral 
measure could not be unambiguously fixed. In a first attempt was made to relate it 
to the gauge fixing procedure necessary to give meaning to the matrix model^, but a clear 
cut derivation of Wyy was not provided. 



In this context the works of refs. 24 1 appeared, where the direct correspondence 



between the Dijkgraaf-Vafa approach and certain generalizations of gauge field theories (on 



non- commutative space-time) was displayed. Immediately after, in ref. |25], the boundary 
condition constraint imposed by the triviality of the AA = 4 superpotential was exploited 
to determine the unknown overall coefficient of the MM functional integral, much in the 
spirit of pq] thus allowing the determination of Wyy. 



similar idea was brought up, in the context of the Eguchi-Kawai one-plaquette model in ref. [ p2[ 
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For the sake of completeness we want to review the derivation of ref . |25[| . There one 
starts with a A^-dimensional MM characterized by the tree-level action (potential) for the 
hermitian N x N matrices , I = 1, 2, 3 

Sm = —tr (g^l [$2, ^-3] + T^aux(<&l) + + ^^i) (2.1) 

and the definition of the Dijkgraaf-Vafa-type free-energy 

iV2 



exp 



9m 



F 

J- rr 



j dl'idl'addaexpf-^^]. (2.2) 



In the above formulae, Qm only plays the role of a scaling constant for the matrix action 
but in the following we will uncover its meaning on the gauge theory side. The matrix 
function VFaux(^i) is an auxiliary potential term chosen according to computational needs. 
For instance, to determine the overall normalization coefficient C^, one has to consider the 
specific form Waux(*^i) = \Mi^\ leading to 

Sm = ^i^( 9^1 [<^2, ^3] + ^<&? + + ^^i) . (2.3) 

9m \ z z z y 

This is nothing else but the so-called MM formulation of the N = 1* model. The A/" = 1* 
model is the = 4 SYM theory supplemented by mass terms for the chiral superfields (all 
belonging to the adjoint representation of the gauge group \}{Nc)). For large values of the 
mass parameters {M/} = {Mi,M2,M3}, this model can be regarded as a regularization 



of the N = 1 SYM theory by means of a mass deformation of A/" = 4 SYM [27, 28 1. 



The N = 1* model is free of UV divergences for arbitrary masses just like the original 



A/" = 4 (superconformal) theory [29, 30, 31, 32, 33 [, as no new UV divergences arise upon 
introducing mass terms. 

The pure N = 1 SYM is reached in the {M/} = Mq — )• oo and g — )• limit keeping 
fixed 

87r2 



A^., = Moexp^-^J, (2.4) 

while in the {M/} — > limit, the AA = 4 SYM theory is recovered for any fixed value of g. 

The key observation in the determination of is that in the {Af/} — t- limit (where 
A/" = 4 SYM is reobtained) one should find that the only holomorphic contribution to the 
superpotential is just given by the tree-level kinematical term. Dijkgraaf and Vafa instruct 
us to compute the effective superpotential by taking the derivative of the planar (which is 
the contribution singled out in the iV — t- oo limit) MM free-energy with respect to S, once 
the latter is introduced in place of the MM coupling constant gm- In formulae, 

W^eff = iVc^^^^%^, F„ = - lim 4log^- (2.5) 

Then, is found requiring that the {-/Vf/} — t- limit of the AA = 1* model in the MM 
setup gives for the free-energy the result 

^-^' = .,lTn^-^'* = ^^ ^ W^eff = 2vrzro5, (2.6) 
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where we have introduced the complexified gauge couphng constant 

47rz 



To = 



9 



2 ^ 2tt 



(2.7) 



Indeed, upon evaluating the matrix integral in ( |2.2D for small {M/}, one gets |25] 



Z: 



Ng 



/ 27rgn, \ 

\NMiM2M^J 



N/2 



[! + ..:_ 



(21 



where the dots stand for terms vanishing in the limit {M/} — )• and = ( ^^^'^ ) ^ is 
the result of the integration over the angular variables in the $-integral^. 

The remarkable fact about eq. ( p.SD is that the ./V-leading contribution to the free- 
energy, F^^^* (see eq. (^^)), does not depend on {M/}. Using eq. ( |2.6[) , one finds 



C 



N 



(2vr)3 63/2^2^ 



7V2/2 



exp 



TTITqN 

Nr. 



(2.9) 



This coefficient, which is a function of N and gm, is assumed to be independent of the 
choice of the potential VFaux(^i) in eq. (|2.1| ). 

Exploiting this result, it has been shown explicitly in 



that choosing Waux(^i) = 
^Mi$2 and in the limit {Mj} = Mq — )• oo, one can also compute the superpotential of 
pure A/" = 1 SYM following the Dijkgraaf-Vafa prescription |17]. The result is 



M/SYM 



-4 



ylog 



3^/2 A3 



1 - log 



A3 



(2.10) 



i.e. precisely the VY effective potential 
2.2 Leigh Strassler deformations 



The Kawai and co-authors formulation of the Dijkgraaf-Vafa correspondence [p3| , 24, P5| 
relies on the finiteness of the = 4 theory. Finiteness is not unique to A/" = 4 SYM, as 
one can concoct various ways (besides adding mass terms ||3^, |3^) in which it is possible 
to deform it without losing finiteness. 

In 1995 Leigh and Strassler Q found that it is not uncommon for a super symmetric 
theory to display a manifold of fixed points, thus discovering new sets of finite theories. In 
particular, they were looking for deformations of the M = 4 model which could mantain the 
beta-function and the anomalous dimensions of the chiral superfields vanishing. Amongst 
them, particular relevance acquired the so-called /3-deformation, which corresponds to a 
modification of the trilinear coupling in the original A/" = 4 superpotential. Explicitly this 
modification amounts to the following replacement 



5tr$i[^>2,^>3] ^ ^tr^>i[^>2,$3]/3, 

[X, Y],3 = e"^/^ XY - e-'^/^ YX, 



(2.11) 



■^We can diagonalize the matrix $ so that J d$ = J Yli dAi ni<7(-^i ~ ^j)^ 



-4- 



where h and f3 are in general complex functions of the coupling constant g. The A/" = 4 
model is recovered at the point P = 0, h = g. 

As we said, the deformed model has attracted a lot of interest after its gravity dual 
was found In particular, it was shown that the /3-deformed model is finite on a whole 
submanifold of the parameter space (/i, /?) 0, 16|. 



If the scalars have a vanishing vacuum expectation value (vev), ($/) = 0, the theory 
is confining and conformally invariant. In the opposite case, for generic values of the 
deformation parameters, it is possible to identify regions of the moduli space (branches, 
in the following) where the gauge group gets spontaneously broken. In order to determine 
these regions we have to solve the F- and D-flatness conditions which here read 

[^1, ^2]/3 = [^2, ^3]/3 = [$3, ^l]^ = (2.12) 

and 

3 

J2[<^i,<^\]=0. (2.13) 
1=1 

Unlike the undeformed = 4 case, in which simultaneous diagonalization of the three 
chiral superfield leads to a solution, now the F-flatness equations are no longer satisfied by 
arbitrary diagonal matrices, but only by special ones. Setting 

{'^i) = dmg{^['\^i'\...,^^il), / = 1,2,3, (2.14) 
it follows that at most only one out of the three numbers cp^'^ , (p^a^ , (p^^ can be non-zero 



for any a S {1, . . . , A'^c}. In this situation it is customary |34] to introduce the three sets of 



integers, Tj, which contain the labels of the non- vanishing vev's. They are 

r, = {a|(^(^) /O}, ae{h...,N,}, 7 = 1,2,3, (2.15) 

and satisfy F/ n Fj = for / / J. Up to gauge transformation, inequivalent branches are 
labelled by the triplet {ni, n2, ns}, where n/ = dim Tj. For the sake of illustration we will 
concentrate in this paper on the particular branches of the type {Nc, 0, 0}, where only one 
of the three chiral superfields (which we take to be $i) develops a vev. 

Although at some special point of the moduli space (corresponding to sets of coincid- 
ing vev's) some subgroup of the original gauge group can remain unbroken, at a generic 



point, i.e. on what we will be calling the "Coulomb branch" [34|, the gauge group will 
be broken spontaneously to U(l)'^''. In this case, the massless spectrum consists of 
"photons", corresponding to the diagonal elements of the gauge field Aa = {V)aa, and 
their gluino superpartners Xaa = {Xa)aa- The two fields combine to form Nc abelian vector 
supermultiplets of A/" = 1 SUSY. We will denote the corresponding field strength by w, 



aa 1 



a = 1, . . . ,Nc. Obviously, there are also Nc massless chiral multiplets, corresponding to 
the fluctuations around the non- vanishing eigenvalues (p^\ a = 1, . . . , Nc- 

What we expect in this situation is that, although no superpotential can be generated, 
the kinetic term for the massless fields may receive quantum corrections. While not much 
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can be said about the kinetic term of the scalars, as the latter is a D-term, the effective 
low-energy action for the gauge fields is a holomorphic function which must be of the form 



W^eff oc ^ r,bu;"''^y;t (2.16) 



a,b=l 



The complex Nc x Nc matrix Tab encodes the effective gauge couplings and vacuum angles. 
Restricting to the {Nc, 0, 0} branch, Tab will obviously depend only on the diagonal part of 
the chiral superfield <I>i, which according to eq. (|2.14|), we shall compactly rewrite as (<J>i). 



At the classical level, Tab is proportional to the identity, t^I = SabTo. Standard non- 
renormalization theorems guarantee that perturbative quantum effects are limited to one- 
loop corrections, while non-perturbative instanton-like terms are expected at any order. In 
other words we will have for Tab{{^i)) an expansion of the type 

oo 

Tabii^l)) = ro6ab + T„V°°''((^1» + ^ rj^^i^l)) e'^''^ = To6ab + fabii'^i)). (2.17) 

k=l 

The one- loop perturbative correction has already been computed in |34, 35 1. 
2.2.1 Different formulations of the Leigh Strassler model 

In the plain = 4 SYM model we have the freedom to (linearly) redefine the chiral 
superfields introducing a pair of conjugate fields. Indeed, starting from the usual form of 
the superpotential tr$i[$2,$3] and defining 

\^ , (2.18) 

the superpotential takes the form ztr <I>i[<I>-|-, <!>_] which up to a z factor has the same 
structure as the one we started from (l.h.s. of eq. ( 2.1l| )) under the formal replacement 
^2,3 — ^ *^+,-- Thus, in the study of the Af = 4 model it makes no difference whether one 
considers the (2, 3)- or the (ib)-formulation of the theory. 

However, when /3 is not zero, it is not true anymore that the form of the superpotential 
is left invariant by the rotation ( |2.18| ). In detail, starting from the (2, 3)-formulation, the 
/^-deformed superpotential in the r.h.s. of eq. ( p. 11 ) transforms under the redefinition ( 2.1^ ) 
as follows 

tr$i[$2,^3]/3 ^tr ($i($^ -$^)sin^/2 + cos /3/2). (2.19) 

In studying the Leigh-Strassler model a choice must thus be made of which one of the 
two possible resulting theories one wishes to study. We can either start from the deformed 
standard superpotential 

^(2.3) =/j$^[$2,$3]^, (2.20) 

which is what we do in the following, or pick up as superpotential 

= (2.21) 
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(obtained by first rotating the fields and then deforming the commutator) which is the 
choice preferred by other authors p6| . As argued above, the two models are different, or 
to say it in another way, rotation ( 2.18| ) and /3-deformation are non-commuting operations. 
Generically, there is no way of morphing one theory into the other. However, we will show 
(see appendix |B|) that, after adding mass terms for the superfields, in particular limits the 
two models become equivalent. 



3. The low-energy superpotential of the massive Leigh— Strassler model 



We now present the results of the MM computation of the effective low-energy superpo- 
tential of the massive Leigh-Strassler model in the confining phase of the theory, i.e. in a 
situation where the chiral field vev's are zero. In the literature, the model has been exactly 



solved in the situation in which the glueball superfield S is integrated out |^, ^ [S^] • 

In the Dijkgraaf-Vafa formulation we are required to compute the free-energy of the 
MM whose action is given by 

3 

. ;\/; r » „ I 

(3.1) 



5™(/3, {Mi}) = — tr J /i<li[l>2, + ^ 

9m I ^ 



While it is clear that, in the limit {M/} = Mq — t- c«, the trilinear coupling becomes 
irrelevant, when keeping a finite (but large) value of Mq we expect to recover the correct 
VY superpotential of pure N = 1 Yang-Mills plus a perturbative expansion in terms of 
inverse powers of Mq. This will guarantee that the Mq — )• oo limit will bring us back to 
the pure M = 1 SYM result. 

The computation proceeds as follows. Starting from 

/Y2 1 c 

d^idisd^se-^™^'^'*-^^), (3.2) 



Z(/3,Mo) = exp 



7? 

2 

9m 



we can immediately integrate over <&2 and <I>3. 
^1, one obtains (see appendix^ for details) 



After diagonalizing the remaining matrix 



Z(/3,Mo) 



(2vr)^g^ 



i i<j 



n[l + 46A2sin2/3/2] ^J]- 



+ e{Xi + Xj-2XiXjCos(3) 



(3.3) 



where e 



lit 



Terms coming from the product over the i index, Y\j^{- . .), in the second line of eq. (^ 
do not yield leading iV^-contributions to the MM free-energy, as shown in appendix We 
are then left with the perturbative expansion of the last product in powers of the small 
parameter e. The dominant term is 



Ko) 



'2T^9m 



2" 



g^9n 



M3 e3/2 



exp 



(3.4) 
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from which the "planar" free-energy immediately follows 



2'^°^M3e3/2 



+ 



Nr 



(3.5) 



At this point, application of the Dijkgraaf-Vafa proposal yields for the low-energy effective 
superpotential the VY result 

o 2 c 



N,S-N,S log 



A3 



(3.6) 



where in rewriting the last equality use was made of the renornalization group running 
of the complexified gauge coupling constant (eq. ( |2.4| )). The result ( |3.6| ) is precisely the 
expected VY superpotential. 

The perturbative expansion of eq. ( |3.3[) in terms of the parameter e gives the higher 
order corrections to the MM free-energy due to finite mass effects. Explicitely, exploiting 
the MM identification of gm with 5, one finds 

Fm-Fl^^ =/\Fm{S,l5,Mo) 

M3 



+ 



) 2 



{b + 2 cos^ /3) 



Ml 



(11 + 12cos^/3) 



3(21-^42 cos2/3-h4cos^^) 



5(^27-M625cos^/3-h440cos^ 



/3) 



(3.7) 



One can check that in the limit /3 — )• 0, i.e. in the pure N = 1* model, eq. (3/7) correctly 
reproduces the known result |^, ^] 



AF„(5,Mo 



Mi' 



+ ... 



(3.8) 



Instead, in the /3 7^ case the contribution to the effective superpotential is obtained by 
adding eq. (3.6) to the contribution coming from eq. ( |3.7| ), and reads (with the superscript 
(2,3) specify that we are deahng with the (2,3)-formulation (eq. ( ^) ) of /3d SYM) 



M/ig''^(5,/3,Mo) = 7Ve5-51og 



A3 



37V, 



h^S 



4c3 



-5^(11 + 12 cos2/3) + ... (3.9) 
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If we eliminate S by the condition which extremizes the superpotential, 



dW, 



(2,3) 



cff 

dS 



0, 



(3.10) 



.(2,3) 



yielding VF^^^g^^^, we can compare our result with the on-shell formula obtained in |^ 
which at all orders in the (zb)-formulation takes the form 



W, 



(±) 



on-shell 



(/3,Mo) 



N,M^ cos 13/2 ei(J3/2\Tn/N,) 

4/i2sin2/3/2 4/i2sin3/3/2 6[{^/2\Tn/N^)'' 



(3.11) 



where 6i{z\t) is one of the standard Jacobi elliptic functions |41| and 



TR = To 



Inh. 



vr 



(3.12) 



As signaled by the superscript (it) attached to the superpotential and as already noted in 



sec. 2.2.1 , the authors of this paper have worked in the (ib)-formulation of the (undeformed) 
= 4 SYM model. Then, in order to be able to make contact with their result, we must 
repeat our previous calculation starting from the expression ( |2.21| ) of the superpotential. 
If we do so, the result for the free-energy differs from eq. (|3.7| ) and reads 



fW = ^5^-Cio, 



Nc ' 2 



g's 

M3 e3/2 



Ml 



■(l-2cos^) + ..., 



(3.13) 



from which we obtain for the effective superpotential in the (zb)-formulation the 1/Mq- 
expansion 



W^j^)(5,/3,Mo) 



N,S-S\og 



9^ 
A3 



Nc 



3iV,^(l-2cos/3) + ... 



(3.14) 



From the discussion above, it should be no surprise that this result is different from the 
one we obtained considering the (2,3)-form of the superpotential, given in eq. (^). 

Now we can finally draw comparisons between the MM computation of the superpo- 
tential and the result of |36|. In fact, starting from eq. ( p. 14 ), we can eliminate S by 
extremizing the superpotential (3.14), finding for the on-shell effective superpotential the 
expansion 



9^ 



5" 



(3.15) 



which coincides with the first few terms in eq. ( 3.11 ). Letting Mq — )• produces a vanishing 
contribution. This is to be expected since in this limit we should recover pure AA = 4 Yang- 
Mills whose only holomorphic contribution to the superpotential is the bare kinematical 
term of the action. 

We have already remarked how, in the /3 — t- limit, the (2, 3)- and the (ib)-formulations 
coincide (and are equivalent to the AA = 1* model). Indeed, if we integrate out S from 
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the effective superpotential (3^) and then send f3 to zero we have the same result we get 
putting /3 — )• in eq. (|3.15| ), namely the Eisenstein series 



nM [q' + 3q' + V + 7q' + 6q'' + 0{q'')] 

gZ 



w. 



(2,3) 



eff 



(/? = 0). 



(3.16) 



Another interesting result (proved in appendix 0) is that, also for non- vanishing /3, in 
the MM formalism the massless limit of the superpotential of the massive Leigh-Strassler 
model does not depend on whether the (it)- or the (2, 3)-formulation is used. 

In conclusion, the web of relations between the two formulations in the various limits 
considered can be schematically represented by the following flow-chart 



(2,3): 



Af-s>oQ 



pure M = 1 




(±): 



/3-j.O 



S out 




/3-deformed 

S out 



■w. 



(23) 



on-shell 



(/3,M) 



M = 1*: 
Eisenstein series 



,9->-0 



■W. 



(±) 



on-shell 



(/3,M) 



4. Spontaneous symmetry breaking in the A/^ = 4 Leigh— Strassler model 

In this section we investigate the behaviour of the massive Leigh-Strassler deformation 
of the = 4 SYM when the gauge group U(A'^c) is spontaneously broken. By varying 
the parameters {/3,Mq) different models are encountered. In particular, when Mq goes to 
zero, we expect to recover the (3d SYM theory, while for /3 — )• we should go back to 
the M = 1* model. Moreover, the limit in which both /3 and Mq vanish must reproduce 
pure A/" = 4 SYM, whose holomorphic contribution to the effective superpotential is, as we 
repeatedly said, just its tree-level kinematical term. As such, it can be used as a sort of 
boundary condition in parameter space useful to constrain the (/3, Mq) dependence of the 
holomorphic superpotential in more general situations. 

If one wishes to induce a spontaneous symmetry-breaking of the form 

n n 

^{Nc)^J{V{Ni), ^Ni = N,, (4.1) 
1=1 i=l 

one has to include in the MM action an auxiliary potential VFaux('^i) of degree n + 1. In 
the N ^ oo limit the dominant configuration in the (zero-dimensional) functional integral 
will be the one in which the matrix eigenvalues are distributed among the n extrema of 
the Waux potential, with Ni eigenvalues located around the i-th extremum. The integers 
Ni satisfy the constraint = A^i -|- • • • -|- A„ and will be also taken to grow infinitely large, 
proportionally to A as A^ — >• c«. 
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For every gauge group factor U(A'j) we have the choice of defining a glueball superfield 
for the full group or for the SU(A^j) subgroup only, namely 



327r' 



327r 



(4.2) 



where Wf is the supersymmetric gauge field strength. Obviously, it is 5^ that describes 
confinement and gaugino condensation, the abelian U(l) degrees of freedom being IR- 
free. However, the form of the Dijkgraaf-Vafa correspondence is particularly simple only 
if expressed in terms of the full glueball superfields, Si, corresponding to the U(A'j) gauge 
groups. In the MM setup, the 5j's are the objects which are identified by the correspondence 



Ni 

Si 44> gm , 
N 



(4.3) 



once the large-iV limit is attained. Thus, for a generic symmetry-breaking pattern, the 
Dijkgraaf-Vafa recipe gives, for the 5-dependent part of the effective superpotential, the 
expression 



i=l 



(4.4) 



while, apart from the tree-level term ToSij, the "coupling constant matrix" of the massless 
U(l)'s degrees of freedom is given by (see eq. ( p. IT] )) 



ds^dSj Ni £^ ^ dSidSk ' 



(4.5) 



k=l 



To obtain the final interesting expression of the coupling constant matrix the last necessary 
step is to extremize the effective superpotential, i.e. solve the equations 



eflf 



dS 



0. 



(4.6) 



The latter is a system of n equations in the n unknowns Si which, once expressed in terms 
of the parameters of the model, will have to be put back into eq. ([4.5|). 

In the following, to simplify formulae, we will concentrate on the branch identified by 
the triplet {Nc,0,0}, i.e. to the case in which the gauge group is U(A^c) broken down to 
U(A^i) X U(A''2) with A^i + = by a non- vanishing This situation corresponds 

to the following particular form of the MM partition function 



exp 



Ar2 



Cj^ I d$id$2d$3e 



-S'™(i>j;/3,Mo) 



(4.7) 



Smi^i; (3, Mo) = — tr [h^i[^2. '^a]/? + ^{^ + ^) + W^^^{^i] 



with 



H^aux(^l) = 7 



^1-2 ('^1 + ^2)$? + 9^1922^1 



(4.8) 
(4.9) 
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where 991^2 ifi 7^ ^2) are the eigenvalues of $1 (see eq. ( 2.14 )) identifying the point on 
the branch we are interested in. The choice of a cubic potential (which has two stationary 
points) allows us to describe the U(iVc) 1— )• U(A'^i) x U(A^2) symmetry breaking pattern. 
Since W^uxi^i) has only the role of inducing the spontaneous breaking of the U(A'^c) gauge 
symmetry, we expect physical quantities be independent of the magnitude of the breaking 
potential. In other words the final answer should not depend on 7. 

After integrating the quadratic $2 and $3 dependence, diagonalization of the remaining 
^1 matrix leads to the formula 



Zn 



-niToN^/Nc 



I' 

h? 

n 



N 



dAj exp 



i=l 

^ + 4A?sin2/3/2' 



J^^aux(Ai 

9m 



-1 



(A. - A,)2 



JJ^ + (A2 + A2-2AiA,cos/3)^ 



(4.10) 



where we have inserted the known expression for given in eq. (|2y 

At this point, the standard procedure consists in expanding each eigenvalue, Aj, around 
either Lpi or ip2 which are the extrema of the potential. We then write 



i = Ni + l,...,Ni + N2 



Ni + N2 = N 



(4.11) 



For convenience, in the following we will set k = i — Ni and pi = qk, whenever i > Ni. 
Then, the symmetry breaking potential becomes 



N 

Waux(Ai 



— -f{fi-f2)Pi-- — 7Pi « = 1,. 

^ 9m "J 9m 

1 ^ / ^ 2 1 ^ 3 , 1 

1^ — l{f2-<fi)qk-o — IQk k = l,. 

^ 9m "J 9m 



(4.12) 



where Ni and N2 are the numbers of eigenvalues chosen to lie around <pi and (p2, respec- 
tively. Defining 



9m 



(4.13) 



and rescaling the pi and qt variables in order to recover the standard gaussian weight by 
putting 



pi^p? = f2 



2 /2 

(Ik Qk 



9m 

7 — {f2 - ^iJQk -Qk, 

9m 



(4.14) 
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one gets 



-niToN^/Nc 



3_ 

h? 



Ni N2 



i=l k=l i<j k<l 



n 



2/3/2] 



N2 



M2 



n + [(9^2 - liqk) e'^'^ -{^2 - lim) e-''^/2][(<^2 - liqk) e"*^/' -(<^2 - ^C^/) e^''/^] 



3((/9i - 992) 



[Y.p'^+^Y.A}- 



(4.15) 



Since in the resulting action the g^'s were displaying an effective negative mass squared, 
we suitably deformed the contour to make the integral convergent, by a sort of Wick rota- 
tion, as is usual in the Dijkgraaf-Vafa approach. Besides, we do not include the contribution 



coming from the product in (the first term of the second line of) eq. ( 4.10 ), as it does not 
lead to relevant terms in the large- ./V limit, as shown in appendix The expression ( 4.15| ) 
will be the starting point of the analysis presented in the following sections. 



4.1 Matrix model perturbative expansion 

The lowest order term is produced by a straight Gaussian integration in eq. ( [4.15| ) , yielding 



^(0) 



-TTlToN^/Nc 



M_^2(gW2_g-./3/2)2 



M_ ^2(^2/3/2 _g-./3/2)2 



2 



M2 + /l2(e«/3/2 _ e-*/3/2 (^2)(e-*/5/2 _ g2/3/2 ^^-^ 
Ni N2 Ni N2 



^1^2 / ^2 

7? 



N4 



Ni 



N2 



n n d^'^ ii(p^ - p^f iK^'^ - ^^^{-\[Zp"+ii a } 

i=l k=l i<j k<l i k 



-TTtTQ 



-«/3/2\2 



Ml 



<^2(g./3/2_g-2/3/2)2 



M2 + /i2(e./3/2 _ e-W2 ^2)(e-./3/2 _ eW2 



N1N2 



h? e3/2 



n1 jv| 

~ / ^,2tCt \~ 



/l2 e3/2 



implying for the leading (tree-level -|- one-loop) contribution to the free-energy the formula 



i7(0) 



1^ 



S? 



(5i + S^f - 5] ^ log 



i=l 



S1S2 log 



h?^{l3)^{-j3)+Ml' 





(4.16) 
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In eq. ( 4.16 ) we have introduced the definitions 



K^ = {-ri 



Ml 



2/^/3/2 



-^^/2^2 



and A{x) = e 



ix/2 



"Pi 



e-"/V2. 



(4.17) 



Following the Dijkgraaf-Vafa prescription we can derive from eq. ( 4.16| ) the coupling con- 
stant matrix for the two massless abelian fields corresponding to the two U(l) subgroups 
of the unbroken U(A'^i) x U(A'^2) g^-uge symmetry. From the formula (see eq. ([4.5|)) [0] 




dSidS2 



^l-loop ^ ^(1) ^ ^ 



(2) 



+ 



(4.18) 



we get from eq. ( 4.16| ) the 1-loop expression 



.1-loop 



log 



M2 + /i2(e*/3/2 _ e-*/3/2 (^2)(e-*/5/2 - e*/3/2 992) 



(4.19) 



which in the Mq — t- limit correctly reproduces the known results of ||3^ , 

The constraint YliNifij = 0, which is automatically satisfied by the definition (^), 
is nothing but the condition ensuring the complete decoupling of the overall diagonal U(l) 
factor contained in the original \J{Nc) gauge group. 

4.2 Next order(s) 



We can carry on our perturbative treatment expanding eq. (4.15) to higher orders in the 
small parameter ^ (eq. (4.13)). One can get in this way the contribution to the coupling 
constants matrix up to n instantons. 

Stopping at second order in the MM formulation (i.e. to order ^^), the complete 
expression of r, upon elimination of Si and ^2 via the standard effective superpotential 
extremization, is given by the quite complicated formula 



^1-loop _^ ^(1) _^ ^(2) 



log 



/i2(eW2 _ e-»/3/2 9j2)(e-*/5/2 _ ^^p/2 (^2) 



+ 



■ [(-3 + 4 cos /3 - cos 2l3)ip\ + 2(-5 + 6 cos /3 - cos 2/3)v9f (^2 



54A4(0)A(/3)A(-/3) 
+ (-13 + 16 cos /3 - 3 cos 2(5)lp\lpI + 8(-1 + 2 cos /3 - cos 2(5)lp\lpI 
+ (-13 + 16 cos /3 - 3cos2/3)(^?v9^ + 2(-5 + 6cos/3 - cos 2/3)99199^ 
+ (-3 + 4 cos /3 - cos 2/3)99^] 
IQh^A^ sin^ 13/2 



+ 



■ [(-139 + 166 cos /3 - 6 cos 2/3 - 26 cos 3/3 + 5 cos 4(3)ipl 



12 



(78A8(0)A2(/3)A2(-/3) 
+ 2(-258 + 368 cos /3 - 125 cos 2/3 + 30 cos 3/3 + 3 cos 4/3)^9^ V2 
- 2(446 - 767 cos /3 + 464 cos 2/3 - 89 cos 3/3 + 36 cos 4/3)^9^ V2 
+ 2(-606 + 1116 cos /3 - 709 cos 2/3 + 250 cos 3/3 + 39 cos 4/3)^9? 
+ (-2513 + 3650 cos /3 - 870 cos 2/3 + 578 cos 3/3 - 125 cos i/3)(pl<fi 
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-4(876 - 1220 cos /? + 871 cos 2/3 + 6 cos 3/3 + 97cos 4/3)v9jy?^ 
+ 2(-827 + 1745 cos /3 - 694 cos 2/3 + 559 cos 3/3 + 99 cos 

+ (vpi o (^2)] (4.20) 
As expected, eq. (4.20) is organized as a power series expansion in the instanton action, 



(4.21) 

and /3 — t- limit, we 
(4.22) 



= exppTTiTo] oc exp[— Svr^/^i^]. 

Remebering that the = 4 SYM theory is recovered in the h? 
can easily check that 

lim T = 0, 

leaving only the overall diagonal contribution, tq (see eq. ( p. 17] )). We stress that, as ex- 
pected, the matrix Tij does not depend on the strength of the symmetry-breaking potential, 
7. This is a quite satisfactory result which confirms our interpretation of the physics de- 
scribed by the formalism as a spontaneous symmetry breaking phenomenon. 



5. Conclusions 

In this paper, we have explored some aspects of the deep connection between matrix models 
and M = 1 supersymmetric gauge field theories. 

We have started considering the massive Leigh-Strassler model, in two of its most 
commonly studied formulations, and computed its low-energy effective superpotential as 
a function of the glueball superfield, for arbitrary (real) values of mass and /3 parameters, 
exploiting the Dijkgraaf-Vafa conjecture. 

We successfully made contact with known results when our formulae are restricted to 
special points of the parameter space, thus providing new evidence for the existence of a 
useful correspondence between supersymmetric gauge theories and matrix models. Along 
the way in this study we have been able to identify the complicated web of relations between 
the different formulations of the Leigh-Strassler model in various limiting situations. 

Then we turned to the study of the spontaneously broken phase of the Leigh-Strassler 
model, introducing an auxiliary potential term to give a non-vanishing vev to one of the 
chiral fields. The MM formalism allowed a detailed and general study of the phenomenon of 
spontaneous breaking of the gauge symmetry. For the sake of simplicity, we displayed the 
expression of the coupling constant matrix governing the dynamics of the left-over massless 
degrees of freedom in the simple case of the \]{Nc) 1— s- U(A^i) x U(A''2) symmetry-breaking 
pattern. We went up to two instantons in the calculation (but there is no problem of 
principles to go higher in the instanton number), finding agreement with the leading order 
(tree- level + 1-loop) expression known in the literature p4| , p5| ]. Reassuringly our final 
formulae do not depend on the "strength" of the auxiliary potential employed to induce 
the phenomenon of gauge symmetry breaking, confirming in this way Dijkgraaf and Vafa 
expectations. 

A challenging open problem which we leave for a future investigation is the extension 
the MM approach to more general deformations of = 4 SYM, such as Tr ^\ and the like. 
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A. Large- iV leading and subleading terms 



In this appendix, we would like to clarify in which sense, in eq. 
we are allowed to forget the term 



+ sin^ /3/2 



), in the large- iV limit 



(A.l) 



Loosely speaking, this is related to the fact that this factor corresponds to a product of only 
N terms, while it is the contribution from the remaining (much more numerous) iV^ — N 
ones that matters. The two terms come from separating the originally unconstrained 
product ■ into an (irrelevant) i = j piece and the set of z 7^ j terms. 

The starting point of the argument is the MM formulation of the massive Leigh- 
Strassler model in its confining phase, whose action in the (2,3)-formulation is given in 
eq. (^T|). The discussion below is limited to the case where Waux is quadratic, but it could 
be extended to a general potential of the form 



W{x) = 



CkX 



(A.2) 



k=2 



Starting from eq. (p.2|), we integrate out two of the three matrices and diagonalize the 
remaining one. After rescaling the eigenvalues by putting 



A,; 



NMq 



(A.3) 



we get 



Z(/3,Mo) = C^J^ 



(2^)^5: 



2^3 

m 



i i<j 



J] [1 + 4eAf sin2 /3/2] ^ H + + " ^^^^^ ^) 



2Ji 



nn 



e = 



(A.4) 
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where the symbols and Yli^j ^^^^ ^^^^ stand for the two factors in the second 

hne of eq. (A.4) and to shorten future formulae we have introduced the "mean value" -like 
notation 

<(*»- / ndA.n(A.-A,f e-iE.A?^. (A.5) 

i i<j 

Our purpose is to compute the logarithm of Z (free-energy) as a power expansion in e, 
according to the formula 



lim ^ log(l + + e^y + 



hm ^ 



eX + e'[Y-^\ + 



(A.6) 



where X and Y come from appropriately collecting terms stemming the various factors in 
the "mean value" (A. 4). 

Dropping an overall constant, which anyway "drops out" when taking the logarithm 
of Z, up to second order in e we find the following contribution to the "mean- value" ( [A.4[ ) 



: n • n » = (( n [1 - ^-^a? + 6eh^xt + oie^)] n 



k,k=ii 



1 - ^ E + E 4 + ^ E E + 0(^^)1 » 



<(i + 4 - Af - 1 a,] + [6.^ + 2/ 5: a? 5: 

i i^j i i ' ' ' 

^ ^ E 4 + ^ E E + E E ^^.] + ^(^') 



k,k^i 



(A.7) 



where we have set Aij = {Xf + A| — 2AiAj cos/3) and introduced the symbol s in place of 
sin/3/2. In the rest of the appendix we will also sometime use C in place of cos/3, so that 
C7 = l-2s2. 

To proceed we need to make use of eigenvalues relabeling invariance in order to make 
use of eigenvalues relabeling invariance to get rid of the various sums in the previous 
equation [42|. Defining 



we can write 



((ni,...,nfc)) = ((Ar...Ar)), 



i 

((^ Aij)) = N{N - l){{Xl + Xl- 2X1X2C)) 
= 2Ar(Ar-l)(((2))-((l,l))C), 



(A.8) 
(A.9) 

(A.IO) 
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so that the coefficient of the first order in e (i.e. X in eq. ( [A.6| )) is given by the expression 

-2s'N{{2)) - N{N - 1) (((2)) - ((1, 1))C) . (A.ll) 
At the next order we have 

{{Y^Xf))=Nm, (A.12) 

i 

i j j 

= N{N- 1){{2,2)), (A.13) 

and 

^ = N{N - 1)^12 J^X- = N{N - l)Au [A21 + (iV - 2)Ai3 

+ (iV - 2)^31 + (iV - 2)^23 + (iV - 2).432 + (iV - 2)(iV - 3)^34] 
= N{N - 1) [Xj + At + 2A?A^ + AXjXlC^ - 4A?A2C - 4AiA^C 

+ 4(A^ - 2){Xi + A?Ai - 2A?A3C + XfX^ + A^A^ - 2AiA^A3C - 2A?A3C 
- 2AiA2AiC + 4A?A2A3C^) + {N - 2){N - 3)(A?A| + A^A^ - 2AiA2AiC 
+ A^A4 + A2A4 — 2A1A2A4C — 2A^A3A4C — 2A2A3A4C + 4Ai A2A3A4C^)j 

= 2N{N - 1){ ((4)) + ((2, 2))(1 + 2C2) - 4((3, 1))C + 2(iV - 2) [((4)) 

+ 3((2, 2)) - 4((3, 1))C + 4((2, 1, 1))C(C - 1)] + (AT _ 2)(iV - 3) [2((2, 2)) 

-4((2,1,1))C + 2((1,1,1,1))C2]}. (A.14) 

In the last expressions, we have introduced the shorthand notation 

N N ff N 

i j i=l j=l ij^j kj^l iy^j k^l 

Lastly we need the formula 

^ ^ A, = NiN - ^ = NiN - l)A,2[\l + A^ + (TV - 2)Xl] 
= N{N - 1) [A:J + A^A^ + {N- 2)XlXl + XjXj + X^ + {N- 2)XlXl 

- 2XIX2C - 2XiXlC - 2(iV - 2)AiA2A2c] 
^ 2N{N - 1) [((4)) + {N- 1)((2, 2)) - 2((3, - {N - 2)((2, 1, 1))C] . (A.16) 



The "mean values" above are, among many others, available in ref. ^2[ where also details 
about their derivation are given. Below we list the ones we need here 

((4)) = l + 2N^ ((3, 1)) = 1 - 2N 

{{2,2)) = 1-N + N^ {{2,1,1)) =2- N 
((1,1,1,1)) =3 

((2))=A^ ((1,1))= -1 
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Using these results, eq. (A^) becomes 



= - lim „ 

N->-oo N"^ 



2 r 

2es2iV2 - eN{N -1){N + C) + e^^(5 + 20^) 
+ e'^N^{8s^ + - 2C^ + 6C - 4) 

+ e27V(4s^ - 4C7s2 - 3C^ + 4C - 1)1 , (A.17) 

where, in order to make clear where each term comes from, use was not made of any 
trigonometric relation. In particular, terms coming from the "diagonal" J^-product are 
all proportional to powers of s^. Looking at eq. ( |A.17 ), and recalling that e is inversely 



proportional to A^, it is possible to appreciate that among the leading terms in there are 
no contributions coming from the expansion of the Y\- term. Moreover, if we now enforce 
the obvious trigonometrical relations, we see that all odd powers of N drop out in the e"^ 
term. The same is true for the even powers in the terms linear in e. 

The contribution to the free-energy can thus be written in the form (see eq. 



Consideration of higher order contributions does not change the conclusion that, in the 
large- limit, we can neglect the contribution coming from expanding the f^-term in the 
expression for the MM partition function. 

B. (ih)- vs. (2, 3)-formulation 

In this appendix we want to show that (the MM formulations of) the massive Leigh- 
Strassler deformations of the J\f = 4 theory, in terms of <^2,3 or ^± are equivalent in both 
the Mo — )• oo and Mq — ;> hmit. 

With respect to the main body of the paper, we will delve in greater detail into 
computations. Our starting points are the following two formulations of the massive /3d 
SYM theories in the MM setup 

z(2'3)(/3,Mo) = /ddidl>2dl'3exp-— tr{/i^i[$2,i'3]/3 (B.l) 

1=2 

zW(/3,Mo) = /"dl>id6+d6_exp-— tr|z/idi[d+,i_];3 (B.2) 

J Sm ^ 

+ Mo6+$_ + iyaux(^l)}. 

We can explicitly verify that these two partition functions indeed correspond to different 
models for generic values of /3 and Mq, by looking at what happens if we perform the 
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substitution 



(±) 



$9 



V2z 



(B.3) 



V2 ' V2 ' 

in eq. ( B.l ). The action, S^'^''^^ in ( |B.1| ) becomes 

5(^'^)(l>i,$±) = — trlMo-l+i)- + «/idi[4>+,6_]cos/3/2 

9m 

+ - 4.^) Sin/3/2 + Praux(^\)}, (B.4) 



showing that generically the two models are different because 



(B.5) 



However, in the hmits Mq — )• and Mq — )• oo, they coincide. To prove this we start by 
computing the partition function of the (2, 3)-formulation. Since the action in (|B.1| ) is 
quadratic in any of the three matrices, one can immediately integrate out one of them, say 
$2- Writing the action in the form [25 



N 



h 



^lA3]}+^^+w^^^{$l)], (B.6) 



2Mo 



we get 



z(2-3)(/3,Mo) = C^ 



27rgm 

NMq 



N r Afn - n 

d$i d$3 exp tr <^ 



(B.l) 



[4>l,|.3]| + W'aux(il)}. 



At this point, as usual, we diagonalize $1, putting (owing to the "gauge" freedom of the 
MM) 

Ai 



$1 = [/ 



Aj G M, [/ unitary. 



(B. 



A 



With the definition (U ^^3U)ij = fiij = fi*^ = [U ^^\U)ij we have 



= -^|AiifcP(A2 + A|-2AiAfcCos/3). 

i,k 

Integration over fiik gives 

Z(2'3)(/3,Mo)=C^J^ 



2vr5„ 



NMq 



l + T72Afsin2/3/2 

_ i_ 

2 



(B.9) 



(B.IO) 



n(^^ - ^^■)' n 1 + ^ + 2A. A, COS /?] e-^ 



-20- 



where has aheady been defined, while the factor ni<j(-^« ~ -^j)^ the usual Jacobian 
associated with the change of variables ( |B.8|) . 

Despite the presence of the M^^ factor in front of the integral, Z^'^''^^ has a finite 
Mq — )• limit, as can be seen by suitably distributing the powers of Mq among the various 
terms. Actually, keeping only the leading contribution (as N — )• oo, see appendix^), in 
the Mo — ;> limit Z^^''^) becomes 



Z(2-3)(;3,0) = C7^J 





"7 


_ Nh . 





i<j 



Xf + Xj - 2X^\j COS ^ 

N 



1 V X — > 

exp{ j;H^aux(Ai)}. (B.ll) 

9m 



If we instead start from eq. ( |B.2| ), we can rewrite S^^\^i, in the form 



Sf\^i,^±) = — tr <j W^aux(^i) + Mo^- 



9v 



Mq 



After integrating out simultaneously the two conjugate matrices <I>±, the resulting partition 
function is 



N 



2-Kgr, 



NMo 



N2 



d$i 



Q 9m 



det 



1 1 + e*/5/2 $1 O 1 - 1 (g) 4>i e-*/3/2 ) 



which, diagonalising can be rewritten as 



zW(/3,Mo) = C^J^ 



iVMo 



n 1 

||l + i|(A.eW 



(Aie^/3/2-Aje-W2) 





"7 


_iVMo. 





/n-iA.n(A.-A,m i^^w,.' 

1 



(eW2 _e-W2) 



n Mo +^(A^ eW2_A^.e-W2) 



Again, redistributing the Mq dependence, we may arrange the last expression in the form 

N2 



zW(/3,Mo) = C^J^ 



2Trgr, 



Nh 



(B.12) 



n 



n - 

.J- + zAi(e*/3/2 - e-W2) ii Mi + .(a^ Q^f}/2 _x. e-»/3/2) 



e Sm 
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Repeating the arguments given in appendix^, and already exploited to deal with the (2,3)- 
formulation, one can show that the diagonal product in the second line does not contribute 
to the large- iV limit. Then, in the Mq — )• limit one gets 



Nh 



-I 



I 'i<j * J ■ J < 



exp< Vl^aux(A,)L (B.13) 

5m ^ 



which coincides with eq. ( B.ll| ), as announced. 

We now move to the Mq — )• oo case. In this limit obviously only the mass terms need to 
be kept in the potential. It is then clear that the two zero-dimensional functional integrals 
are just gaussian and identical. 

Actually, the latter is quite an expected result, since, in the dual gauge theory descrip- 
tion, it is known that in the Mq — )• oo limit we decouple completely the ^2,3 degrees of 
freedom in one case and the <I>± degrees of freedom in the other. What is left is simply an 
integral over the last chiral field that is not affected by the rotation. In other words, there 
is no difference between the two possible forms of Z(f3, Mq) in the Mq — )■ 00 limit. 
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